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Abstract 

We consider AdSs and warped AdSs vacua in new massive gravity and study the highest weight modes 
and general propagating modes as a set of solutions for the linearized equations of motion. We observed 
that depending on the choice of gauge there are two types of solution. We show that for warped AdSa 
vacuum, the massless modes which only appear in harmonic gauge have zero energy density and they do not 
get higher curvature corrections. By computing the energy density it can be shown that all massive modes 
have negative energy density. Our computations prove that the massless modes for squashed space-time are 
always stable, while for the stretched space-time there are restrictions on the parameters of the theory to 
have a stable mode. 
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1 Introduction 



According to AdS/CFT correspondence [T] the boundary fields parametrize boundary conditions of tlie bulk 
fields and couple to the operators in the dual conformal field theory. The sub-leading radial behavior of these 
fields at boundary is obtained by finding the most general asymptotic solutions to the field equations. For 
theories that admit asymptotically locally AdS solutions these general solutions, which is sometimes called 
the Fcfferman- Graham expansion, can always be found by solving algebraic equations [5]. Even though higher 
derivative terms are treated as perturbative corrections to two derivative actions, they do not change the usual 
AdS/CFT setup [3]. In fact the field equation implies that there are some additional boundary data to be 
specified. 

Since the three-dimensional gravity described with Einstein-Hilbert (EH) action has no degrees of freedom 
[4]l5], commonly other deformations of pure EH gravity have been considered. In fact the addition of higher order 
derivative terms provides the theory with propagating degrees of freedom, i.e. three-dimensional gravitons. The 
first studied theory was the topologically massive gravity (TMG) which constructed by adding a cosmological 
constant and a gravitational Chern-Simons term |6] - [8] . The other theory which we are going to consider here, is 
the New Massive Gravity (NMG) [S| . Also some extended theories of NMG have been discussed in [TU] , [H] , [H] . 
The quantization of such theories seems to give a richer structure than the EH theory and provides interesting 
toy models for higher-dimensional theories of quantum gravity. 

The AdSs vacuum is the first known solution of three-dimensional massive gravity models. Another vacuum 
solution of the higher order derivative actions is the warped AdS^. These vacua admit black hole solutions 
known as BTZ [8] and warped AdSs black holes [131 [14]. In this paper we are going to discuss about the 
asymptotic behavior of the metric fluctuations in NMG model. Since the corresponding black holes have 
equivalent asymptotic behavior as vacuum solutions we consider only the vacua here. 

To find the behavior of metric fluctuations around a background there are two approaches. In the first 
approach one finds the highest weight modes which are belong to some representations of the isometry group 
of the background. According to the isometry group of each background {SL{2, R)l x SL{2, i?)^ for AdS^ and 
SL{2,R) X U{1) for warped AdS^) there are some generators. The highest weight modes are those which are 
killed by raising ladder operators (Li and Li for AdS^ and Li for warped). Additionally these modes must 
be the solutions of the linearized equations of motion. Of course there is no reason that all solutions to the 
equations of motion fall in the highest weight representation. 

In the other approach, one substitutes the metric fiuctuations as eigen-modes of energy and momentum in 
the linearized equations of motion and tries to find a decoupled differential equation for each component of the 
perturbations. The regular solution at the asymptotic limit can be obtain by using the Frobenius's method. 
In the asymptotic limit where the radial direction, r, approaches to the boundaries, solutions behave like , 
where the parameter i? is a function of constants in the theory and it can be shown that it is closely related 
to the frequencies in the highest weight approach. These asymptotic solutions often called general propagating 
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modes. In the case of TMG the calculations for warped AdSs using these two approaches has been discussed 
in [IS]. 

The stability around a certain background depends on the selection of consistent boundary conditions. For 
example there exist several consistent choices of boundary conditions for three dimensional massive gravity 
models in AdS^ background. But only more confining ones exclude negative energy modes (unstable modes), 
which arc different to Brown - Henncaux (BH) boundary conditions |16| . Similarly the massive propagating 
modes of warped AdS^ do not obey the Compere - Detournay (CD) boundary conditions [TTlfTS] . In fact the BH 
and CD boundary conditions are only consistent with the pure large gauge in which the asymptotic perturbations 
are given by the Lie derivative of the background metric h^^ ~ 'C^g^^ and ^'^ is the non- vanishing asymptotic 
diffeomorphism. We have been considered these boundary conditions in [11] which lead to the central charges 
of the dual CFT living at the boundary of asymptotically AdS^ space-times. Note that for the propagating 
modes the metric perturbations cannot be written in this form for any ^'^ of background metric g^i,. 

The stability of TMG around AdS^ vacuum in harmonic gauge has been shown in [T^] . They have obtained 
a stable theory without negative energy at chiral point fil = 1. Calculations for stability of BTZ black hole in 
NMG has been done in [21], [22] . 

In this paper we will study the perturbations around AdSs and warped- AdS's backgrounds in NMG for two 
different gauge conditions, the harmonic (transverse) gauge given by h'^'^ = and the other gauge which is 
given by /i^ip ~ and check the stability of different modes. 

This paper is organized as follows: In section 2 we bricfiy discuss the Lagrangian of NMG and its equations 
of motion and extract the linearized form of these equations around an arbitrary background. In section 3 
we review and study AdS^ vacuum and its perturbations in two mentioned gauges, then check the stability of 
this vacuum. In section 4 we repeat all steps in section three but for warped-iris's vacuum. In section 5 we 
discuss the stability conditions for warped solution in different gauges. We compute the energy density for the 
corresponding fiuctuations. In Section 6 we will consider extended NMG model which contains sixth order of 
derivative terms and study its effects on different modes, section 7 includes summary and discussions. 

2 NMG and its linearized equations of motion 

The new massive gravity is described by the following Lagrangian [9] 



where A is the cosmological constant. In addition to a three dimensional gravitational coupling constant, there 
are two other coupling constants, ki = and K2 = which are specified by a mass parameter m. The 

equations of motion corresponding to variation of gravitational field are 




(2.1) 



T, 



nNMG 



= R^u - ^gtJ.uR + ^9,.u + 2KiR{Rf,^ - -g^,^R) + (2ki + K2){g^,^a - V^V^)i? 



(2.2) 
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Since our goal in this paper is the study of behavior in the gravitational fluctuations, we need to linearize the 
equations of motion around an arbitrary background. Let gf^^, be a background metric and its fluctuation is 
described through g^^, = g^^ + h^^. The linearized Christoffel connection, Riemann and Ricci tensors are 

1 



per 
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y^Vt^K + VaV,/i;j - y^V.h - V„V"/i^,) , 6R = g"^ SR^p - h"^ R, 



(2.3) 



All derivatives are taken with respect to the background metric 5^1/ . Using these relations the linearized 
equations of motion for NMG Lagrangian become 

^j.WA/G ^ gj^^^ _ ^h^^^ji _ 2A) - ig^, SR + 2Ki(i?^, SR + RSR^,, - ^ V i?^ - i g^,R5R) 



— 2hP^ Rfj^p^cr R-l — -rhpuRta + 7,9fiuh^^ R\Rp<j — -cji^^RP'^dR, 



" -•- ^ pa- 



pa) 



(2.4) 



Everywhere we use the bar notation, it means that a quantity must be computed by background metric. 

3 AdSs vacuum in NMG 



The global AdS^ metric is a vacuum solution in new massive gravity 



ds^ = V 



-{l+r')dT' + 



r^dip^ 



(3.1) 



1 + r2 

where I is the radius of AdS^ space. The isometry group of this space is SL{2,R)]^ x SL{2,R)ji with the 
following left and right moving sets of generators 



\/l + 



(3.2) 



Each set creates a SL{2, R) Algebra, = 2Lo , [Lq, L±i] = T-^±i- 

Inserting AdS^ metric (|3.ip as a background into the NMG equations of motion (|2.2p . fixes the value of 
cosmological constant in terms of the mass parameter and radius of AdS^ space 

4m2/2 + 1 



A 



(3.3) 



4to2^4 ■ 

In order to study the linearized equations of motion and find a solution to these equations in a specific back- 
ground, we need to fix the gauge freedoms. In what follows we will consider two different gauge fixing conditions 
and compare their results. 
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3.1 Harmonic gauge 

The AdS^ metric due to its symmetries makes the hnearized equations of motion (|2.4p more simpler. For this 
metric the curvature tensors are 

1 2 6 

R^pva — J2 id tJ,i^ 9 prr 9p<y9pv) j Rpu — j29lJ.f J R= -j^ ■ (3-4) 

If we compute the trace of hnearized equations of motion by multiplying equation (j2.4p by g^'^ and impose the 
harmonic gauge condition V^/i,,^ = then we will find a traceless condition for fluctuations, h ~ /i^^ ~ 0. 
The harmonic gauge and traceless condition simplify linearized equations of motion ()2.4p to a fourth order 
differential equation as follow [19] 

(^(i)cW^(M)^(M)/^)^^ = . (3.5) 
The covariant derivative D defines as a first order operator 

= (5/ ± ;e,"^V„) , (I?(M/M))^^ ^ ± ie^^^/^v J , (3.6) 
where e^f^'' = "^^^^ = Lagrangian one finds fi as PU]. ^T] 



. x/2 + 4m2/^ 

^ = 21 ■ (^-^^ 

The above differential equation describes one left and one right moving massless graviton mode in AdS^. There 
are also two massive degrees of freedom. These can be seen by the following equations 

(i?(^)7^(^)/.)^. = P(V^-^)V=0, 

p(A/)^(M);,)^^ = ij(V2 - I + M')h,, = , M' = p:'-^ = ^^^^^^ . (3.8) 

The mass square relation shows that in order to have stability (Tachyon free condition), the parameters of the 
theory must be bounded to rn^P > i values [21] . 

3.1.1 Highest weight solutions 

It is possible to find solutions for linearized equations of motion by requesting the highest weight conditions. In 
other words these solutions belong to the representations of isometry group of AdS^ background. The highest 
weight solutions are eigcn-modes of SL{2,R)l x SL{2,K)ii generators with the following eigen- values 

L+ih^n = , L^h^i, = — ^ — h^jj , Lj^ih^y ~ , L^h^i, = — ^ — h^i, . (3-9) 

This leads us to consider the following ansatz for metric fluctuations 

-(l + r2)/i(r) 51 (r) g^ir) 

V = e-*("*-'^^)g(r) I g^{r) {l + r^r^f2{r) gz{r) |. (3.10) 

92{r) gsir) r'^ fsir) 
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The conditions L_|-i/i^jy = L+i/i^^ = hold when fc = 0. By solving differential equations from the highest 
weight conditions, the following behaviors for unknown functions in the ansatz can be achieved 

fi{r) = -^{Ci + Csr^ + Csr*) + K^Ci + C2 , 
1 + z 

r(l + r^j 2 

93ir) - -JT^r^AC^r^ + {iC, - HC^ + 2C:i)r^ + iC^ - ZiC^ - C4 + 2C3) , (3.11) 
r(l + r^) 

where we have used the traceless condition to fix /3(r) = — /i(r) — f2{'f')- In this solution, Ci's (z = 1, . . . , 5) 
are constants of integrations and function g{r) must be 

g(r) = (1 + r2)-^r'= . (3.12) 

The harmonic gauge /i^i/ = fixes three of C^'s, as 

^ _ »(fc + l)(fc-2c^ + 4) ^ ^ _ »(fc^ + (3-c^)fc-c^^ + 2c^) ^ i^k + u^ 

2(k-uj + 2){k-^ + i) 2{k-Lo + 2)ik-L, + 3) ^''^'-AikTYr'- ^^"^^^ 

Gathering all these results in the linearized equations of motion gives six equations for two unknown constants 
Ci and C3. One can easily verify that for every arbitrary subset of equations there is a relation of the form 

^2x2 ) =0. (3.14) 

To have a non-trivial solution (Ci 7^ 0, C3 ^ 0) the determinant of A4 must be zero which gives a relation 
between the frequencies co and wave numbers k 

w = fc,fc-2,fc + 4,fc- 1 ± i V2 + 4m2/2 , fc + 3 ± i V2 + 4m2/2 . (3.15) 

But there are some constraints which must be considered here. First of all as we mentioned earlier, the highest 
weight condition holds for fc = 0. On the other hand we expect that the asymptotic metric fluctuations falls off 
faster than the background metric. We also demand the Tachyon free condition, m'^P > i. 

Therefor if Ci 7^ then w > 4 and we have two possible solutions: Either w = — 1 + ^ V2 + Am?P for 



TO^p > or w = 3 + iV2 + 4m2Z2 for mH'^ > \. 

If Ci = then all constants except C3 are zero and we need w > 0. Here we have four possible solutions: 
w = 4,3 + i\/2 + 4m2Z2 and w = -I + ^V2 + Am'^P for m'^P > i and a; = 3 - ^\/2 + Am'^P for ^ > m'^P > \. 

3.1.2 Asymptotic behavior of propagating solutions 

To find the asymptotic behavior of metric perturbations let us to consider the following ansatz 

(l+r2)/i(r) 51 (r) g^ir) 

/i^, = ;2e-M-M I g^^r) -{I + r^Y' f2{r) 53(r) |. (3.16) 



Since we need to know the asymptotic behavior if we drop or keep 1 in 1 + r'^ the final rcsuhs as we have 
checked will be exactly equal. So for simplicity of calculations one may drop 1. In harmonic gauge, we have 
three relations among the unknown functions in /i^^. These are 

r^g[ir) + 3r^giir) + iujr^ fi{r) + ikg2{r) = , 

-rVsW - 2r/2(r) + ^fc.93(r) + iLogi[r) + r{h[r) + h[r)) = , 

r^g'^{r) + ir^ gz[r) - ikr^ h{r) + iug2{r) = , (3.17) 

so we can write fiir), /^{r) and g2ir) in terms of gi{r),gs{r) and /2(?') and their derivatives. On the other 
hand we have g'^'^ST^J^^'-^ ~ which together with the harmonic gauge condition implies traceless condition for 
metric fluctuations. Therefore we have /i(r) + /2(^) + /3('') = 0. This constraint fixes another function and 
one finds 

53(0 = ^('■'/2(0 + 3r/2(r) - jc^gi(r)) . (3.18) 

If in equations of motion we impose the gauge condition p.l7p and constraint p.lSp then we will find six 
differential equations for two unknown functions gi{r) and f2ir). Among these equations the rr component is 
a differential equation for /2(?') 

2rVi^^ + SGrVa" + ((191 - 2m^fy - 4(fc2 - cj^y)/;^ + {{-UmH^ + 329^ + 20{uj^ - k^y)/;, 
+2((-8m2z2 + my + {mH^ - ^^){k^ - u?y + (fc^ - u?f)j2 = . (3.19) 

In order to analysis this equation and find its asymptotic behavior we use the Frobenius's method. Let's insert 
a series solution f2{r) = r^^ S^^o above equation and then take the r — > oo limit and insert the 

coefficient of the greatest power of r to zero. In this way the possible values for B are (cq ^ 0) 

B = 2,4,3±iV2 + 4m2Z2. (3.20) 

What about the other unknown function gi{r)l We observed that one cannot write a linear combination 
of equations of motions in order to find a differential equation for gi{r). We need to consider higher order 
derivatives of equations of motion as well. After doing some calculations we will find an 8th order differential 
equation for gi{r). Since it is a very lengthy equation we have not write it and just present the final results. 
The asymptotic behavior of this differential equation can be found again by assuming gi{r) = r^^ X^J^o 7^- 
The values of B' are 

B'=:l,3,3,5,2±iv/2 + Am^P ,A±^V2 + 4m^P . (3.21) 

It is possible to find a relation between /2(f) and derivatives of gi(r). By this relation one can verify, which 
values of B' are related to values of B. We find that for every value of B there are two values for B' in such a 
way that 

B' = B±l. (3.22) 
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Again wc demand that the asymptotic metric fluctuations fall off faster than the background metric. This will 
be happened if B and B' > 0. 

Consider B' = B + 1 then when B > automatically B' > 0. In this case all solutions of B are acceptable 



without any constraint except one. The only constraint happens when B = 3 — i-\/2 + 4to^P which is acceptable 
for i < m2/2 < IT., 

When B' = B — 1 then B > 1. Again all solutions are acceptable and the only constrain is coming from 



B = 3- iV2 + 4TO2p ^here i < m^P < |. 
3.2 hf,^ = gauge 

We are interested to know the behavior of the theory under another gauge fixing condition. So we do all steps 
in previous section but now for a new gauge h^^^ = 0. In this case it seems that the equations of motion cannot 
be written in the decoupled form (|3.5p . Despite this, we can study the solutions similar to previous sections. 

3.2.1 Highest weight solutions 

The metric perturbations ansatz are given in p.lOp . The highest weight equation gives the following functions 
(note that we have not traceless condition in this gauge) 



JiKr) ^ -[C2 + iCi)r — — h Ci 

1 + 



(1 + r^j ^ ' 
h[r) = r C2+ iCiT + iiCi + C2 + Ci H j- , gi{r) - 



J.2 ' ?'(1 + r^) 

r(l + r2) 



g2{r)^{C2+iCiy + {2iCi + C2y-iC^+iC^. g^{r) ^ . (3.23) 



To impose the gauge condition /i^^ = in the linearized equations of motions we have to choose Ci = C2 = 
C4 = C5 = Ce = 0. Now there are six equations for just one unknown constant C3. The only consistent 
solution to these equations is C3 — 0. Such a solution is called the pure gauge solution. In other word, to have a 
highest weight solution with nontrivial values for frequency and wave number in this new gauge, we must have 
fiir) = Qiif) = for i = 1, 2, 3. The same behavior previously has been observed in [5T] for BTZ black holes. 

3.2.2 Asymptotic behavior of propagating solutions 

Although we showed that the highest weight solution in the new gauge /i^j^ = is a pure gauge, but let us try 
to find the asymptotic behavior of metric perturbations in this gauge. We start from the following ansatz 

/(l + r2)/i(r) .91 (r) 0\ 

V = ^'e-*^"*"''''M 9i{r) H^ + r^yf2{r) . (3.24) 

\ y 

To find an asymptotic solution, we insert the above ansatz into the equations of motion (to simplify the 
calculations we have ignored 1 in 1 + , the final results are exactly the same). There are five independent 
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equations of motion which arc not enough to decompose the differential equations. On the other hand we have 
g^^ 5T^J^^'~^ = which imphes 5R = 0. If we use this equation together with its first and second derivatives then 
we can obtain a differential equation for fi{r). If insert ml — ^ then for /i(r) we find the following differential 
equation 



+ r-^)^(^+^^)r^ + 3(e^-^f(fc^---2)r* + 3(^+?^)(fc^-o.^)V)fcVi=0. (3.25) 

If we consider an asymptotic series solution for the above equation, i.e. fi{r) = X]J^°=o 7^' then we can read 
the values of B from equation B{B^ — 6B + ^ — C^)(5^ ^ 25 + \ — = 0. So the possible values of B are 

B = 0,B=l± i + 2 , S = 3 ± ^^4^2 + 2 . (3.26) 

Doing the same computation but for /2(f) = r^^ 7^ gives another complicated differential equation. 

The final result for values of B' will be 

-B' = 2,l±iV4C^ + 2. (3.27) 

For gi{r) = r^^ J2'^=o 7^ ^^"^ ^^"^ following values of B" 

B" = l,B" = 2± ^^4^2 + 2 , B" = ±^ ^4^^ + 2 . (3.28) 

To suppress the blow up of metric fluctuations in asymptotic region, we must limit power behaviors to -B > 0, 
B' >0 and B" > 0. 



4 Warped AdS^ vacuum in NMG 



Another vacuum solution to the equations of motion with less symmetries than the AdS^ background is the 
warped AdS^ vacuum 



ds^ 



We define the warp factor as 



(l + r2)dr2 + — ^ + ^—(rf^ + rrfr)2 . (4.1) 
I + t/^ + 3 J 



2j^ , ^ 

(4.2) 



The generator of time translations is given by a Killing vector as ( ^'i^^ [~1 + •^{^71^)1"^])^^ dr- For i/^ < 1 it 
is always time-like but for > 1 there is a transition surface from time-like to space-like. For < < 1 the 
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warp factor < cr < 1 and the space-time is called squashed. For > 1 wc have 1 < a < 2 and it is called 
stretched. The special case u = a ^ 1 corresponds to a deformed AdSs with a fibration [T5] . 

The isometry group of AdS^, i.e. SL{2,R)l x SL{2, R)ii, is broken to SL{2,R)l x C/(1)_r due to presence 
of the warp factor a. The generators of this symmetry transformation can be constructed out of the Killing 
vectors and are given by 

r „ . r T„ 1 



Ln 



idr , L 



±1 



±e 



-dr =F iVl + r^dr 



VT 



J 



-id, 



v ■ 



(4.3) 



If we consider the NMG Lagrangian and insert the vacuum solution (j4.ip into equations of motion p.2p 
then we will find the following values for NMG mass parameter m and cosmological constant A in terms of the 
warp factor a and the warped AdSs scalar curvature {R = —-p) 



3 21cr^-4 



A = - 



3 21(7^ _ 72(T^ 



16 



(4.4) 



2^2 cr2 - 4 ' 2/2 21cr4 - 88cr2 + 16 ' 

Since the warped factor is limited between < ct < 2 we can specify the sign of parameters of the theory for 
different values of the warped factor. Supposing P > (or negative curvature) then we have the following 
behaviors for mass parameter and cosmological constant 



a 


< 0- < (Tl 


CTl < (7 < (72 


(72 < (7 < (73 


(73 < (7 < 2 
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< 


> 


> 


> 


A 


< 


> 


< 


> 



Table 1: Behavior of mass parameter and cosmological constant in NMG for < cr < 2. 

The critical values in warped factor which have been appeared in the above table are 

2^189-42^1! 



0-1 



0.436 , (72 



0.489, (73 = ^^l^^±i^^ 1.786. 



(4.5) 



/21 21 21 

In what follows we will find different behaviors of metric fluctuations around the warped AdSs metric for 
different gauge conditions. 

4.1 Harmonic Gauge 

As we mentioned before since the warped AdSs background has less symmetries than AdS^, we do not have the 
simple rules as equation p.4p . This makes the equations more complicated, so in most cases we just present 
the final or important results. 

4.1.1 Highest weight solutions 

Unlike the AdS^ vacuum, here we consider the behavior of metric fluctuation only in presence of the harmonic 
gauge condition V^/i^'' = and we have not a traceless condition. The highest weight condition for metric 
perturbations are given by the following relations 



J h^y = k h^v , L^h^v = ojh^y , Lihfj,v = , 



(4.6) 
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where L±i,Lo,J are Killing vectors that generate the SL{2,R)fi x U{1)l isometry group which are given in 

(|4.3p . Like the previous case we choose the following ansatz 

/ Mr) /2(r) /3(r) \ 

/i^.(r,r,^) = /6(rK('=^--) Mr) h{r) h{r) . (4.7) 

V h{r) h{r) Ce / 

The (3, 3) component of the highest weight condition Lih^^, = fixes the value of fe{r) to 

/6(r) = e'=*-""'-(l + r2)-^. (4.8) 



Using this value, and by solving differential equations from other components of Lihfn, = one finds the 
following solutions 

/i(r) = Ce - C3 - 2iC2 r + C4 + iCi r - C5 , /2(r = , 

2(1 + r"^) 

Cs r^ + C^r + C3 Cir+_C2 
(l + r2)2 ' •'^^'■^ 1 + r 

In the case of harmonic gauge we can fix three of the Q's in terms of other constants, i.e 



/3M-(C6-»C2)r + zCi, /4(r) = , hir) = ^\2X ■ (4-9) 



■ 2 

= ^^^|^-jy(C4fc2 + 2((a2-l)C3 + (l-c.)C5)fc + a2C4(2-c.)), 



^ 2( 



^6 = 4^ ( - 2^3^' + ((t' + 2c^ - 3)C74fc2 + ( ((-2C5 + 2C3)c^ + 2C5 - 4C3)a2 



^^^l^-jy (2C3fc2 - C4((t2 - (2 - c^))fc + 2(C3 - (1 - w)C5)a2) 
1 

2fc(fc2~f cr4) 

+ 2(l-cj)((w-l)C5 + C3))fc-a2C4(l-c^)(2-cj))(T^ (4.10) 

We now substitute the ansatz (|4.7p in the linearized equations of motion and use the above gauge conditions. 
To have a non-trivial solution, the determinant of coefficients of C3, C4 and C5 for each subset of equations must 
be zero. Independent of choice of subsets we find always two polynomials P\ and ¥2- Defining w ~ + ^ 

Pi = — (2(16A:^ - 136fc2 + 9)cr^ + (64fc^ - 368fc^ + 652^^ _ 45)cr^ - 4A:2(48A:^ - IGO/c^ + 95)cr'' 
64 ^ 

+ 16fc*(12fc2 - I9)cr2 _ 64fc6-) ^ J_(4(4fc2 _ 5^8 _|_ (-43^4 _ -^52^2 + 59)cr^ - 8k^(12k^ - 17)cr^ 

' 16 

+ 48fcV) u + + ^(4fc2 - 5)a'^ - 3fcV4)u2 + ^6^3 ^ q ^ 

P2 = 84cr*^ + 24(2fc2 - 5)(7^ - — (176A:4 + 424^^ - 333)cr^ + i(176A:^ - 208fc2 + 9)a^ + {Uk* + -)k^ 

16 8 2 

+ (a2(48a4 + ^cr2- 5) -2/c2(licr4-12cr2 + i))u_ (72(11^2 _ 2)^2=0. (4.11) 
Solving these equations gives the frequencies of allowed modes. We will discuss these modes in the next sections. 
4.1.2 Asymptotic behavior of propagating solutions 

In order to find the behavior of metric perturbations at the asymptotic limit as a propagating mode let us 
consider the following ansatz 

l2/A „2\ 



/2/4__ v^ , -(l + r2)/i(r) cr2gi(r) ra^g2ir) 

r<J^92{r) cr2g3(r) -0-^3(0 
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(4.12) 



If we impose the Harmonic gauge we will find three equations between unknown functions as 

(1 + r2)2/^-i(l + ((fc(a2 - l)r^ + a^ior - k)g^ - a\kr + u)g^ - ir{h - h)) -ra\{r^ - 1)52 - h) =0 , 

0-2(1 + r'^fg[ - i{l + r2)(fcr + - r(^i{{a'^ - VjUr"^ + a'^ujr - k)g2 - 2a^{l + r^)gi^ , 

a^{l + r'^fg'^ + i{{kr + w)rcr2 - fc(l + r'^))!^ + a'^r{i{kr + uj)g2 + 2(r2 + 1)33) = . (4.13) 

If we find fi{r), /^{r) and g2if) from above gauge conditions and insert them into the equations of motion we 
will find six mixed differential equations for three remaining functions. To find the asymptotic r — >■ 00 behavior 
of the solutions there are two approaches. In first approach similar to previous cases, we recombine equations 
of motion and try to find a decoupled differential equation for each unknown function. In the second approach, 
we can use the following behaviors for the remaining functions 

gi(r)^Cir-^, /2(r) ^ Car"^ , g3(r) ^ Csr"^-^ (4.14) 

By putting these relations into the equations of motion and by going to the asymptotic region one finds six 
equations for three unknown constants Ci, C2 and C3. To find _B, it is enough to select three out of six equations 
and insert the determinant of the coefficients to zero. These values of B again must be independent of the choice 
of equations. Doing all steps above, we will find again two polynomials from highest weight approach (|4.1ip 
just by replacing w by B. 

4.2 h^^ = Gauge 

The calculations in this section is roughly analogous to one accomplished for TMG in [TS] while we do it here 
for NMG. 

4.2.1 Highest weight solutions 

The situation here is very similar to the previous case. Wc start from highest weight conditions in (|4.6|) and 
find exactly the same value for fG{r) as ()4.8|) . Substituting it in the remaining equations we will obtain other 
functions as 

3C2r2 + 2{Ce - C^)r - C2 



Mr) ^ (Ci - Cey + 3C2r - C3 , h{r) 



2i(l + r2) 



hir)=Cir + C2, h{r) = (^^^2)2 ' - ,(1 + ^2) ■ (4-15) 

Next, we use the gauge fixing condition ft,^^ = such that the highest weight condition is preserved. We choose 
C'e = /al^) ~ /sl^) = which states that Ci ^ C2 — 0. Wc now insert this perturbations into the linearized 
equations of motion. We can show that for each subset of equations of motion we have a matrix 

(4.16) 
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where the propagating modes are those with detM — while the pure gauge modes obtain by detA4 7^ or 
equivalently C3 = C4 = C5 = 0. There is a unique factor in aU determinants of M for aU subset of equations 
and we find that detA4 ^ P2. So in this gauge just one of the previous polynomials survives. This polynomial 
gives the value of possible frequencies for each propagating mode in the highest weight 

- + ^ [96 cr^ + (31 - Uk^) - (10 - 48A:2)cr2 - 4fc2 ± a{ 1500cr^^ 
2 V11't2-2L V ; V ; ^ 

+ (409 + l%k^)a^ - (68 + AOk'^)a^ + (4 + 26fc2 + k^)a^ - 2(2 + k^)k^a^ + k^^ 

96(7^ + (31 - AAk^) - (10 - \'$,k^)a^ - ik^ ± 4f ISOOcr^^ - 1116cr^ 



Wl,2 
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1 1 

W3,4 = : 



2 VTT^ 

+ (409 + 18/c2)ct» - (68 + 40fc")CT« + (4 + 26fc^ + fc*)tT^ - 2(2 + k'^jk^a' + ' . (4.17) 



4.2.2 Asymptotic behavior of propagating solutions 

To find the behavior of metric perturbations in h^^ ~ gauge let's consider the following ansatz 

/ -(l + r2)/i(r) Mr) 0\ 

V = e-'("*-'='^) /2(r) (l+r2)-i/3(r) . (4.18) 

\ / 

As before, we have six equations of motion where just five of them are independent. One the other hand we 
can use the trace of energy momentum tensor g^^ 5T ^j^i, = 0. If we use this equation together with its first and 
second derivatives we will have eight equations. After a little computation we can find decoupled differential 
equations for /i(r),/2(r) or fsi^r). For each function if we write an asymptotic expansion series, we will find 
the same behavior. The fall off power again is given by a polynomial, B = + 5 

P3 = 84(7^ + 24(2fc2 - 5)a^ - — (176A:^ + i2ik^ - 333)a* + -{176k-^ - 2Q8k^ + 9)a^ + (life'' + -)k^ 

16 8 2 

+ (cr2(48cr4 + ^cr2- 5) -2/c2(llcr4-12cr2 + l))u- (72(11^2 -2)^2=0. (4.19) 

As we see P3 is equal to P2 in the previous gauge (j4.1ip . As another equivalent approach to find the asymptotic 
r — > 00 of the solutions, we can insert the following behaviors for fiuctuation functions 

h{r)~^Cir-^, /2(r)^C2r-^, ^r) C^r^^ . (4.20) 

By inserting these values into the equations of motion and by going to the large values of r one finds six Algebraic 
equations for three unknown constants Ci,C2 and C3. The only consistent nontrivial solution happens when 
we have the P3 polynomial for values of B. 

5 Stability 

There are two main conditions that must be checked in order to have a stable solution. The first one is the 
positivity of the energy of a typical solution and the second one is the reality condition for frequencies. In this 
section we perform both these checks to find the domain of reliability of our solutions. 
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5.1 Energy condition 

To find the energy condition we follow the same approach presented in [TJ] for TMG (one may also use the 
ADT construction ^^). According to [23] the conserved charges Q{£,) of a diffcomorphism invariant theory, 
which are associated with the Killing vectors are given by the following expression 

Qio - / HeEi^M'^]dxn , (5.1) 



16nG 7e 

where E is a spatial hyper-surface at constant time and * represents the Hodge star operation. To find the 
conserved energies, we substitute the first order perturbations h^/}} of the highest weight solutions, into the e'^^J , 
the energy-momentum pseudo-tensor. The energy density of the gravitational wave is given by |15j 

^~/^^^/-Ir^4l^^, (5.2) 

where we have used ^'^ ~ (1, 0, 0) to find the energy. We consider the physical perturbations in their real form 

-0,11/ = a hi_i,^ + a* {hi_,y)* , (5.3) 

where /i^i^'s are given by (|4.7p - (|4.9p in the harmonic gauge and by (|4.7p together with (|4.15p in the other gauge. 
After applying the gauge conditions we can remove the coefficients Ci, C2, Cg and write C3, C4 in terms of C5 
from linearized equations of motion, i.e., 

6*1(0;, fc,(T) <j'^ 5*2(0^, fc,g) 

^3 = '"?"? — I — t'-^S' ^4 = -77 7r7 — r — r <-5 , (5.4) 
iio(a;,/e, cr) k bo[uj,k,a) 

where So, Si and S2 are real functions and arc given in the appendix A. 

To avoid the divergences we consider the energy density per unit length in the ip - direction as in |15j . The 
final result for energy can be written as the following sum 

e = \aC^\^Y.\\ n , 2W+4 [b^XKo)) ^ Mk,a)BUk,a). (5.5) 

The above integrals are finite for Re{ijj) > i and for n < 8 and they obey a recursion relation [T5] 

(2w + 7 - n)A„ = 2fcA„_i + {n - 1)^.-2 , (5.6) 

this enables us to write the energy in terms of Aq which is real and positive. As we showed in previous sections 
there are two types of solutions labeled by two polynomials: 

• Massive modes: 

If we insert the P2 polynomial which exists in both /i^,^ ~ and harmonic gauges we will find a negative energy 
density in < cr < 2 and for all values of k. This negativity of energy for massive modes is similar to the 
TMG [T5]. Therefore all massive modes in the squashed or stretched space-times are unstable. 
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• Massless modes: 

Interestingly if we examine the energy density of modes corresponding to the Pi polynomial which exists only 
in the harmonic gauge, the result will be zero. So this means that these modes are describing the massless 
modes of the theory. For TMG in [TS] since they have considered only the h^j,^ = gauge they have not seen 
such massless modes. To find these modes in TMG we preformed computations of [15] in the harmonic gauge. 
Again the results contains two type of polynomials. The massive modes found in the /i^j^ = gauge and another 
mode. The new mode has exactly the same polynomial structure as Pi in NMG case and its energy density is 
zero. 

5.2 Frequency condition 

As we mentioned before another check for stability of solutions is the reality condition for frequencies. In order 
to do this we need to solve the polynomials and find their roots. These roots were the values of B (fall off 
powers for propagating solutions) or they were the allowed frequencies lu for highest weight modes. We demand 
that these roots to be positive and real valued. 

Since we used a change of variable as a; = ±it5 + i and because the energy condition restrict us to Re{uj) > i 
so only the plus sign is allowed and the reality condition for frequencies translates to Reiu) > 0. 

Although the P2 polynomial shows a negative energy modes but in what follows we will check the regions 
of real frequency for massive modes as well. 

5.2.1 P2 polynomial (Massive modes) 

Starting from P2, we find the following roots for u 

"± = A^2,t^t'\. . ^1 = ^'(96^' + 31^' - 10) - 4(a^ - l){lla' - l)^ , 
4cr^(llcr^ — 2) 

A2 = ct'^(4 - 68a^ + 409a'^ - 1116cr*^ + 1500tT^) + 2a^{9a^ ~ 2){o^ - \fk^ + (ct^ - \fk'^ . (5.7) 

As we mentioned before, to have a real value for w, we must have positive real values for u±. According to the 
roots of A2 we have two situations: 

• When A2 has two real roots for fc^ 

1. To have a real valued uj^ we need A2 > 0. Suppose that solving A2 = in terms of fc^ gives two real roots 
so we can write A2 = (cr^ — l)2(fc2 — I]_)(fc2 — S+). So the reality condition of u±_ will restrict fc^ to (note that 
E_ > E+) 

(T^ — 1 

2. On the other hand A2 itself must be real, which means that Ej- must be real. This will be possible if we 
choose < cr < w in the admissible interval of < cr < 2. 
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3. One must notiee that in the interval < fi < y ^ = 0.426 we have always S]_ > 0, but for 0.381 < u < 
0.398 we have < 0. So in this sub-interval we have just > S_. 

4. By looking to the values of u±, in the interval < a < \ we see that the denominator has a negative 
value for both u±. So u± > if Ai ± 4A| < 0. In this interval Ai changes sign from negative values to 

positive values and it is a monotonically increasing function for each values of k'^, and — 4fc^ < Ai. Also 4A| 

i_ 

is a positive monotonically decreasing function for each values of k^ and 4A| < 4k'^. 

i_ 

5. Finally we observe that |4A| | > |Ai| for all values in the mentioned interval and for all fc^, so we conclude 
that just w_ is an acceptable solution. 

• When A2 has not real roots for fc^ 



1. In this case we are restricted to y ^ < tr < 2 where A2 has not any real roots and we have A2 > 0. So for 
it± > we need Ai ± 4A| > 0. 

2. In Ai = a'^{96<7^ + Slcr^ - 10) - 4(cr2 _ l){na^ - the first term (coefiieient of k'^) is always positive 
in < a < 2 but the coefficient of k^ changes sign from positive values to negative values at tr = 1, so Ai 

can be either positive or negative. 

1 

3. If I All < 4A| or equivalently Aj - I6A2 < then is allowed but u_ is not. On the other hand if 

1 

I All > 4A| or equivalently Ai — I6A2 > then we have two choices. If Ai > then both u± are allowed and 
if Ai < then neither m+ nor u_ are not allowed. 

4. Consider Ai - I6A2 > 0. If < ct < 1 then always Ai > so 

' 4(.^-l)(ll.^-l) ' ^'-'^ 

but in this interval Sg < so fc^ > and therefore there is no restriction here. But 

Al - I6A2 = a^iUa^ - 2) ( - 3cr^(448CT^ - 640a'^ + llltr^ + 6) - 8{a^ - 1)(96(t'* + 430-^ - 9)k^ + I76{a^ - i fk'^) , 

(5.10) 

which by Af — I6A2 > assumption k^ restricts to 

fc2 < S+ , fc2 > S_ , S± = —4 -(96(T^+43cr2-9± (81 -576o-2+3784cr4+24000cr**- 128640-6) 5) . (5.11) 

44((T^ — 1) 

In this interval we find that S_|_ < Sq < S_ so S_|_ < therefore fc^ < S+ is not allowed. If we look at S_ we 
see that it changes sign from positive values to negative values &i a — 0.506. The final result is as follows: 
T 



< cr < 0.506 ^ S_ > ^ fc2 > S_ , 0.506 < (7 < 1 ^ S_ < ^ /c^ > 0, (5.12) 

and two solutions u± are allowed. 

5. Consider Ai — I6A2 > 0. If 1 < cr < 2 and Ai > then Sg > fc^. In this interval of a we have 
always S_ < Sq < Here Sq > and 2+ > but S_ changes sign from positive values to negative values 
at (7 = 1.103. So two solutions u± are allowed, either when k"^ > S+ in all the interval or if k'^ < S_ when 
1 < CT < 1.103. 
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6. Consider Af — I6A2 < where just u+ was allowed. For ^ < a < \ wc find ^+ < < ^_ but since 
in this interval S+ < and S_ changes sign we conclude that for < a < 0.506 we have < /c^ < S_. 

For 1 < (T < 2 we find S_ < < S+. Here 2+ > and S_ changes sign so for 1 < cr < 1.103 we have 
S_ < A:^ < S+ and for 1.103 < cr < 2 we have < A:^ < S+. 

5.2.2 Pi polynomial (Massless modes) 

The Pi polynomial was au^ + bu^ + cu + d = where 

a = , b = ^4(3(^2 - l)fc2 + ^^2(2^2 _ ^5)) ^ 

c = ^2(3(^2 - Ifk* + ^a^a^ - l){2a'- - I7)k^ - ^^'(20<72 - 59)) , 

d^{a^~ l)3fc6 + 1(^2 - 1)2(2^2 - 19)(72fc4 _ ^(74(^2 - l)(68a2 - 95)fc2 + ^^7*^(2^2 - 5) . (5.13) 

As we see from these coefficients for < cr < 1 we have always a > 0, b < 0, c > and d < 0. According to 
Descartes' rule of signs this polynomial has either three or one real positive root. So at least there is one real 
positive root in < cr < 1 for all values of k^. 

Dividing Pi by and changing variable asw = t+ | — ^ — " ^^ k^ we can write the polynomial in its 
depressed form, + pt + q = Q 



p a'{a' + l2) + 12{a^-l)k' -cT^(a^ - 36) + 18(^2 ^ l)(17a^ „ 6)fc2 _ 

^ 12cr2 108cr2 ■ ^ ' ' 



The roots of this equation depend on the sign of its discriminant A = ^ + |^. IfA>0 then there is just one 
real root, t = (—f A^)^ + (— | — A^)'^ and if A < then there are three real roots for t. The discriminant is 

(0-4-4)2 (cr2-l)(3cr«-cr6-98cr4 + 36cr2 + 24),, (863cr4-612cr2 + 60)(cr2-l)2 (0-2-1)3 

A = /CH -. k — k . 

432 216ct2 4320-4 270*^ 

(5.15) 

It is very hard to find exactly where we have a change sign in A but since cr is limited between and 2 we can 
compare values of k with cr in three regions. 

• k^O 

In this case only the first term in A is dominant (A = — '''^432 ) so for all values of cr its value is negative and 
we have three roots for t. In this region u = t + j + ^ and t has three positive roots t ^ i ' I ' t ~ 

• few 0(cr) 



In this region there are sub-regions where we have positive or negative values of A. For example if we consider 
k = a then A = j - + ^o-^ - 5a-^ + y^a^ - j^a'^" where we have two roots at o- = 0.60 and a = 1.05 
between these roots A < and beyond them it is positive. 



16 



• fc > cr 

At this limit the last term in A is dominant and we have A ^ ^ ^'^27a^ ^° ^® have a changing sign around 
the cr = 1. For cr < 1 there is one real positive root and for ct > 1 there are three negative real roots. All roots 
are near the value of u = — fc^^— S-i. 

For clarifying the behavior of roots for u we have presented some numerical values in the following table. 



cr 


1 

3 


2 
3 


4 
3 


5 
3 


fc = 


0.25,2.25,1.19 


0.25,2.25,1.03 


0.25,2.25,0.36 


0.25,-0.14,2.25 


fc = 1 


9.97 


2.00 


2.26 


2.30 


fc = 2 


33.9 


5.90 


1.67 


1.39 


fc = 3 


73.9 


12.1 


0.23 


-0.83 


fc = 10 


802 


126 


-34.6 


-52.8 



Table 2: Real roots in different region of parameter space. 



These numerical values show that for < cr < 1 we always have at least one positive root where by increasing 
fc only one positive root remains. But for 1 < cr < 2 by increasing fc the number of positive roots reduces and 
for large enough values of fc there is no positive root in this region. 

6 Extended new massive gravity 

In this section we consider higher curvature corrections to NMG and try to find their effects on spectrum of 
massless and massive perturbative solutions for both previous background metrics. The Lagrangian of extended 
NMG up to third order curvature terms [TO]- [12] is given by 

^ENMG ^ ^(^^3 ^R'^PR^^ + KiRR^.^R'^'' + K5 R^) , (6.1) 

where K3 = — 3^ , K4 = and K5 = — . The equations of motion are obtained by the following energy- 
momentum tensor 

+ K4(i?M.i?a/ji?"^ + 2RR^''Ra,-^g^,RRc.pR''^ + +,g^..V„V^(i?"''i?) - 2V„V(,,(i?,)"i?) 

- [V^V, - .g^,n](i?„^i?"^)) + K5(3i?p.i?' + 3[.9^,n - V^V„]i?2 - . (6.2) 

In what follows we will study and solve the linearized form of (|6.2p around the AdS^ and warped vacuum. Since 
the behavior of the extended NMG in different gauges is exactly similar to NMG we present only the results in 
the harmonic gauge. 
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^ = SZnre • (6-3) 



6.1 Ad S3 vacuum 

Similar to NMG wc can find different properties of AdS^ vacuum solution in Extended NMG. In this case the 
equations of motion restrict the cosmological constant to 

imH'^ + 2m^P + 1 

Once again if we perturb around the AdS^ vacuum we will find a similar fourth order differential equation (j3.5p 
but in this case 

— ^^^?pti — ■ ^^-'^ 

The differential equation describes again a massless graviton mode and a massive graviton with mass square 

so the Tachyon free condition occurs for m^P > or — ^ < rr?!^ < . 

The highest weight solutions can be found by inserting the ansatz p.lOp into the highest weight equations 
p.9p and determining h^^i, . Substituting these values into the linearized equations of motion accompanied with 
the gauge conditions gives the following values for frequencies 



where the highest weight conditions restrict us to fc = 0. Comparing these results with p.isp shows that only 
the mass of the massive modes gets correction. 

We can confirm the above result by looking to the asymptotic behavior of the propagating solutions by 
inserting the ansatz p.l6p into the linearized equations of motion and using the gauge conditions. As an 
example we obtain the following differential equation for /2 (r) 



4 4-1/ 

2 



_ ((7^4;4 _ ^^2;2 _ 651 ^ ^o^^2;2 ^ i)(^2 _ ^2))^3^/ ^ ((_8^4;4 ^ g8,^2p ^ 33)^4 

2 8 2 

+ {mH^ ^m^f |)(fc2 - + {m'P + l)ik' u:^) h ■ (6.7) 

The series solution around the boundary behaves as such that B has the values 



1 8m^l^ + 47772/2 I 3 

Again the relation between uj in the highest weight mode and B in the propagating mode is very similar to the 
NMG case p.20p and only the square-root term has been corrected. 

6.2 Warped- AdS's vacuum 

If we consider the NMG Lagrangian and its curvature corrections and insert the vacuum solution (|4.ip into 
the equations of motion T^^^*^ + T^/^^*^ = 0, then we will find the following values for mass parameter and 
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cosmological constant (A = y^Ula^ + 264cr2 - 80) 

3 210-2 -4± A 



2 2 

m = m. 



A = A+ = 



4/2 0-2-4 

-7830-s + 23400-* + lOOSo^ - 576 ± (-GSo" + 2160^ - 48) A 



(6.9) 



(o2 -4)(21o2 -4± A)2p 

By taking ^2 > Q the following behaviors for mass parameter and cosmological constant will be obtained. The 
reality of cosmological constant constraints the warp factor to begin from Oc = ^V-627+342\/6 ^ 0.509 . 



a 


Oc < cr < 


OL < 0- < 0_R 


an < C7 <2 


m\ 


> 


> 


> 


A+ 


> 


< 


> 



Table 3: Behavior of m+ and A+ for o^ < c < 2 . 



ol and (Jfj are two real (between and 2) roots of A+ = 0. Numerically they are equal to o^ = 0.558 and 
(Tji = 1.802. 








715 <^< 73 


73<-<2 


9 

m_ 


> 


< 


> 


A_ 


> 


> 


< 



Table 4: Behavior of m_ and A_ for Oc < o < 2 . 



6.2.1 Massless and massive modes 

Similar to NMG we use the ansatz ()4.7p and put it into the equations (|4.6p . We achieve the same results as (|4.9p 
and (|4.10p . If we substitute these results into the linearized equations of motion then we would find non-trivial 
solutions when the determinant of coefficients is zero. Like the NMG in harmonic gauge, here we have two 
polynomials P[ and P'2- The first polynomial is P[ = Pi, in other words we find again the massless mode of 
(|4.1ip and it does not receive any correction. The other polynomial which represents the massive modes gets 
correction and it can be written as 

P2 = Eo{k,cr) + Ei{k,<T)u + E2{k,a)u^ + E3{k,(T)u^ = 0, (6.10) 

where lj ~ ui + ^ and the functions Ei are given in the appendix B. 

If we try to find the asymptotic behavior of propagating solutions again we will find the above results exactly. 

7 Summary and Conclusions 

In this paper we have mainly discussed about the stability of AdS and warped AdS vacua as solutions of the 
new massive gravity in three dimensions. First we found the equations of motion for NMG and then linearized 
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them around arbitrary background 5^^. We determined the value of cosmological constant for each kind of 
solutions in terms of other parameters in the theory such as m,l,a. 

We observed that in the harmonic gauge one could write the linearized equations of motion as a product of 
decoupled operators p.Sp for the AdS vacuum. We should note that although there have not been found such 
a decomposition for the warped space our calculations might be a hint for finding it. 

We considered the behavior of metric perturbations in the asymptotic limit from two points of view. In the 
first view, fluctuations satisfied in the highest weight equation in addition to the linearized equations of motion. 
While in the later, these equations decomposed to some differential equations for each perturbation that could 
be solved by analytical methods in the asymptotic limit. We have used two different gauge conditions to write 
the linearized equations of motion. 

In warped AdS^ vacuum the value of frequencies uj, coming from the highest weight formalism, was related 
to the radial fall off power parameter of the general propagating modes at the boundary. In fact in most cases we 
have studied, the propagating modes of metric fluctuations were also belong to the representations of isometry 
group of the background. 

We observed that the existence of a mode in a vacuum, depends on the gauge choice. For example in the 
case of warped vacuum in the harmonic gauge we obtained two polynomials, Pi and P2, while in hf^ip ~ gauge 
there was just one polynomial P2- So it seems that the behavior of propagating solutions or highest weight 
modes depends on choosing of the gauge conditions and some modes are forbidden when we change the gauge. 
This can be seen for AdS^ vacuum as well, since as we found in harmonic gauge there is a highest weight mode 
while it becomes pure gauge when we go to the h^^ip ~ gauge. 

We also discussed comprehensively about the stability of these vacuum perturbations at the asymptotic 
limit and found the domains of validity for parameters in different gauges. We showed that Pi polynomial 
describes the massless modes of the theory. For the squashed warped space-time there is always one possible 
allowed frequency but for the stretched warped space-time the stability is limited to some regions of space of 
parameters. We also showed that P2 which describes the massive modes has always negative energy both in the 
squashed and in the stretched warped space-time. 

We have considered the extension of NMG, which are constructed from curvature terms and are consistent 
with the AdS / C FT context. This extended Lagrangian did not change the main results but only corrected the 
values of mass for the massive propagating modes of AdS^ background. 

In warped space-time we showed that again there arc two polynomials. A polynomial P{ in the extended 
theory which is exactly equal to Pi so that the massless modes do not correct by higher curvature terms. In fact 
this encouraged us to look for this result in the TMG model. Again we found two polynomials in the harmonic 
gauge which one of them was exactly Pi and the other one was the result of [H] for TMG in /i^^ = gauge. 
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A 5*0, ^i, 5*2 functions 

• The real functions in h^^p = gauge 

50 = 42(a;-l)CT^" + (15w^-45w2 + 3(5A:2 + l)w + 48 + 6fc2)CT*^-(2w^ + 7a;'* + (4fc2-54)w3 + (9/fc2 + 68)w 
+ (2/s'*-44/s2-8)cj + 50/s2+2fc'*)CT^ + (2cj^-10w^ + (4fc2 + l6)w4_(8fc2 + 8)w3 + (2fc4 + i7fc2)^ 

51 = 42(w-l)a^V(15w^ + 27w + 15/fc^a;-51a;2 + 24)CT^-(2w^-10tj'^+(4fc2 + 31)a;^-(25fc^ + 68V2 
+ (2/c^+4U-2 + 60)w-18A:^-15/s**)cr^ + (2a;5-(2fc2 + 10)u;'*+(16fc2 + 16w3-(8 + 4fc'*+40A:^)w2 

+ {51k^ + Uk^)uj^36k^^23k'^~2k^)a^ + 2k\4+k^-'iLu+Lu^){{2k^+uj^-2Lo)a^-k^), 

52 = 42(fc^-w(a;-2))CT^-(15a;''-45w^-6w^ + (57fc^ + 72)w-36fc^-15fc'^)cr*' + (2a;^-8a;^ + (6 + 2fc^)a;'^ 
+ {8-26k^)uj^ + {45k^-2k^-8)Lj^ + {'i0k^-18k^)-2k^ + 19k'^-60k^)a^uj+4:{uj^-5uj^ + {7+2k^)uj^ 
+ (^-4fc2)cj2 + (Zfc2_7^^4)^_^^4_2fc2)fc2CT2-4(a; + i)(4+fc2-4a;+w2)fc4, 
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• The real functions in harmonic gauge 

= (42fc^-(42w2-420a; + 336)A:2-84a;^ + 84w-294w2+294w3)cr" + (78fc^ + (99w2+417w-420)fc'^ 

-(36w'*-624a;^ + 630w^ + 996w-1056)fc2-57w^+606a;2-117a;''-375w^ + 207tJ^-264w)cr^2-(10fc* 

+ (28w2-200a; + 296)fc6-(24w'*-398a;3 + 542a;2 + 1142w-1034)fc''-(4a;6-196a;^ + 266a;'^+900a;3 

-1368w^-548cj + 1104)fc2+2w**-2w^-20w^-2w5 + 160w^- 134^3- 196w2 + 192w)cri° + (4fc^" + (16w2 

-32w + 30)fc* + (24cj''-100w''' + 168w^-462w + 398)fc'^+(16w^-108w^ + 270w'^-664w^+497cj2 + 1021cj 

-912)fc'* + (4w^-44w'^ + 156w^-242a;^ + 56a;'*+450w3-606w2-56w + 384)fc2-4a;2(w2-l)2(a;-2)3)cr^ 

-(16fc** + (52w2_ii2tj+46)fc^+(60w''-264cj3 + 348w2-324w + 156)fc'^+(28w^-192w5+426w''-318cj3 

-154cj2+408w-192)fc2+4w(w^-l)(w^-6w^+4w + 6)(w-2)2)fcV^ + (24fc^ + (60w^-144cj + 58)fc'^ 

+ (48w''-228a;3 + 304w2-62w-88)fc2+4(3Lj'*-9w'''-3w2 + llw+4)(w-2)2)fc''a'' + 4(fc2 + (Lj-2)2)fc8 

-4(fc2 + (w-2)2)(4fc2 + 3tJ^-4(w + l))fcV, 

S*! = (84fc''-168fc2-42a;'^ + 84w3 + (42fc2+42)a;2 + (210fc2-84)w)ai''-(78w6-333a;^ + (99fc^ + 522)a;'^ 

-(876fc2 + 537)w3-(36fc'*-1443fc2-534)cj2-(543fc'^ + 600fc2 + 264)w-57fc*'+828fc'* + 240fc2)CT^2 

+ (10w^-58a;'^ + (28fc2 + 80)w^ + (28fc2 + 110)w'5 + (24fc'*-196fc2~292)w'^ + (230fc'^-22fc2-22)w3 

+ (4fc^-702fc'* + 270fc^ + 364)w^ + (144fc^-376fc'^-588fc2-192)a;-2fc^-426fc^ + 860fc''+792fc^)cr^° 

-(4cji°-28w^ + (16fc2 + 64)w®-(84fc2 + 24)w^ + (24fc'^ + 154fc2-l08)w^-(84fc'^ + 114fc2- 132)^5 

+ (16/c^ + 132fc''-14fc2 + 8)w''-(28fc6-144fc'^-294fc2 + 80)w^ + (4fc^ + 58fc^-518fc'^-474fc2 + 32)a;2 

+ (234fc^ + 191fc'^-72fc2)w + 16fc®-585fc^-124fc^+384fc2)cr** + (16cj®-100w^+4(13fc2+48)w^ 

-4(57fc2 + 5)w^ + (60fc'* + 290fc2-304)a;'*-(156fc'^-6fc2 + 216)w^ + (28fc^ + 136fc'*-194fc2+96)w2 

-4(7fc^-32fc'*-38fc2 + 24)cj + (4fc^ + 38fc''-232fc2-176)fc2)fcV^-(24a;^-132w5+4(15fc2 + 52)cj'' 

-4(51fc^-5)w^ + (48fc'' + 178fc^-256)w2-(72fc''-54fc2 + 80)w + 12fc^+36fc''-80fc2+64)fcV'' 
+4(3A:^+4a;2-3w-4)(fc2+4-4w+w2)fc6cr2_4;,8^^2_^4_4^_^^2^^ 

52 = (252w3-252w2 + (252fc2-252)w-504fc2)(7i2^(270w5-330w'* + (540fc2-702)w3 + (426-1080fc2)w2 
+ (744+270fc'*-516fc2)a;+984fc2-750fc'^)cr^"-(24a;'^-262w6 + (712 + 72fc2)w5 + (270 + 382fc2)w''-(72fc^ 
+ 1564fc2_868)a;^-(22fc'^-1434fc2 + 340)w2-(24fc^ + 852fc^-992fc2+480)w + 1652fc''-456fc^-142fc^)CT^ 
+ (8w^-52w^ + (108 + 32fc2)a;'^-(152fc2+40)a;^ + (48fc'* + 216fc2-116)w'5-(144fc'^ + 340fc^-96 
-108fc'* + 892fc2 + 36)a;3-(40fc^-40fc'^ + 356fc2+40)w2 + (8fc**+954fc'*-668fc2)w+4fc* + 340fc^-1086fc^ 
-24fc2)a6-(24w^-132w^ + 8(9fc2 + 26)tJ^-4(63fc2 + l)cj'' + 8(9fc^ + 19fc2-22)w3-2(54fc'*-31fc2- 12)^2 

+ (24yfc^-56fc'' + 156fc2+48)w + 12fc^ + 270fc'^-216fc2)fcV'' + (24w^-108w'^+4(12fc2 + 29)tj3-4(24fc2 
-ll)w2+8(3fc'^-2fc2-7)cj + 12fc'^+88fc2-24)fcV2-4(2cj + l)(fc2+4-4cj+cj2)fc^. 
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B Functions of P2 polynomial 

£;o=[-(430272CT^2+(232704fc2-84096)cr^"-(70608fc'^-285672fc^+696405)CT^-(2048fc^+176160fc'*-573856fc2 
+ 158790)a^+(6144fc^ + 129488fc'^-6392fc2-32544)CT'*-(6144fc^ + 12928fc'* + 65216fc2-2880)CT2-1100 
+ 16128fc2 + 2048fc^)A + 5614272cri'*+(3034368fc2 + 3341952)cri2_(923472/(;4_587e236fc2 + 13288653)cr^° 
+ ( 30720fc^ + 1409376fc'*-5291328fc2- 1399494) CT^- (116736fc^-1164336fc''+7862728fc^-3083592) cr^ 
+ ( 165888fc^-3181376A:'* + 5537120fc2-381312) - (104448fc^-1849600/c^ + 1422592^2 + 18432) cr^ 
-129024fc2_3l8464fc4+24576A:^]/5614272, 

Si=[(-29088CT^°+(17652fc2~56073)cr^-(768fc'^ + 15280fc2 + 5218)cr^ + (2048fc'^-8484fc2 + 8000)cr'' + (7584fc2 
-1792fc'*-1744)cr2_i472fc2 + 512fc'*)A + 379296(T^^-(230868fc2-1034721)cri°+(11520fc'*-3072fc2 

+410238)cr^ + (478852/c2-255560-35840fc'^)cr^ + (43264fc'^+328944A;2 + 5488)cr'^ + (98112A;2 + 10816 

- 25088fc'')cr^ - 14080fc^ + 6144fc'^] /701784 , 
£;2= [(4413(T* + (3370-384fc2)cr6 + (896fc2_i952)cr4 + (-640fc2 + 640)(T2 + i28fc2)A-57717(Ti"- (89622 

- 5760fc2)cr* - (13952fc2_ 27272)cr^ + (1472 + 12160fc2)cr'^ + (-3328- 5504fc2)cr2 + i536fc2] /350892 , 
£;3 = 32(150-2+4- A)cr2(cr2 - 1)2/87723, 
where A = ^Jllla'^ + 264 cr2 - 80. 
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